Modeling and simulation of flotation reagents system in anionic reverse iron oxide flotation at different temperatures

Removal of quartz from iron ore was accomplished industrially via an anionic reverse flotation technique. However, in that kind of flotation, the interaction of the flotation reagents with the components of the feed sample makes the flotation a complicated system. Thus, the selection and optimization of regent dosages at various temperatures were performed using a uniform experimental design to estimate the optimum separation efficiency. Besides, the produced data as well as the reagent system were mathematically modeled at different flotation temperatures, and the graphical user interface GUI of MATLAB was conducted. The advantage of this procedure is that the user interface displayed in real-time can be conducted by adjusting the temperature at different values to automatically control the reagent system, besides predicting the concentrate yield, total iron grade, and total iron recovery.

Iron ores are the primary source of iron which is essential for the world's iron and steel industries and serve as the backbone for the country's basic infrastructure improvement. The iron ores consist of iron oxides mainly of magnetite (Fe 3 O 4 ) and hematite (Fe 2 O 3 ). However, the average grade of these ores in China is less than the cut-off grade of 45% 1 . Thus, upgrading process development was required for increasing the grade of the ores to meet the demand of the iron and steel industries. The beneficiation was studied extensively 2-4 which consists of several stages, and anionic reverse flotation was the common route in the processing plants to produce the final concentrate [5][6][7][8][9] . On the other hand, direct flotation as an alternative route led to a partial flotation of quartz and in consequence reduced the grade of the final concentrate 10 . However, the anionic reverse flotation is established at a pH value of 11.5 to maximize the repulsive electrostatic force between negatively charged iron oxide (valuable) and the quartz (gangue) particles. Also, a depressant such as corn starch is utilized to inhibit the floatability of iron oxide particles besides aggregating fine iron oxide particles to not be transported to the froth zone 11 . In addition, an activator such as calcium oxide is used to activate selectively the quartz by changing its surface charge into positive. Then an anionic collector like commercial TD-II is adsorbed on the positively charged quartz to increase its hydrophobicity and in consequence its floatability [12][13][14] . Thus, the appropriate reagent system in anionic reverse flotation is very critical and complicated especially when dealing with low-grade iron ores as a result of the interaction of the different flotation reagents with the diverse components of the feed slurry. Thus, the selection and optimization of the reagent system are vital investigation stage which needs a lot of effort and is time-consuming, thus a uniform experimental design has been conducted by the authors 1 . Some of the advantages of uniform test design are a very quick, highly efficient, economical experimental design method that reduces test times, shortens test cycles, and is capable of quickly finding multi-factor optimization schemes 15,16 .
The anionic reverse flotation of iron oxide is usually carried out with automatic slurry temperature control at 35 °C 12,14 or 30 °C 1 . In this paper, the flotation experiments of the iron oxide feed were carried out using the uniform test design method at different flotation temperatures from 20 °C up to 45 °C to investigate and optimize the influences of the flotation reagents as a function of flotation temperature.
In this type of research, various diagrams can be quickly generated using a MATLAB-based graphical user interface GUI. It is one of the fast and most common software used for data analysis and visualization, algorithm calculation, and numerical design 17 .

Methods
Materials. The representative experimental samples were collected from the flotation circuit feed of the Anqian processing plant located in Anshan, Liaoning. The industrial plant runs a flotation circuit with samples having 90%, 50%, and 10% finer than 104 μm, 35 μm, and 6 μm respectively; and total average iron content of about 48%. The complete sample characterizations were illustrated in previous work using X-ray diffraction (XRD), X-ray fluorescence, and Mineral Liberation Analyzer (MLA) 1,12,14 . The major minerals are iron oxide (hematite/magnetite) and quartz with other impurities such as pyrite and apatite 1,12,14 .
The flotation circuit is conducted based on the typical plant conditions using the plant chemical reagents including starch as iron oxide depressant, lime (CaO) as quartz activator, and commercial anionic TD-II as quartz collector. Besides, analytical grade NaOH and HCL (15% w/w) were used as pH modifiers.
Flotation process. The studies were conducted using the plant flowsheet shown in Fig. 1 but at different flotation temperatures from 20 °C up to 45 °C. The closed circuit consists of roughing, cleaning, and three scavenging flotation stages at 11.5 pH value 1 . At the roughing stage, depressant, activator, and collector are added, while www.nature.com/scientificreports/ a third of the collector is added in the cleaning stage. The tailing of the cleaner is combined with the concentrate of the first scavenger and sent back to the rougher feed, the concentrate of the second scavenger is recycled back to the first scavenger feed, and the concentrate of the third scavenger is sent back to the second scavenger feed (Fig. 1). The products of the flotation circuit were dewatered, dried, homogenized, weighed, and the total iron grades were analyzed.

Uniform test design.
For the flotation circuit of iron oxide, a uniform test design was conducted to estimate the optimum flotation reagent systems at different flotation temperatures from 20 °C up to 45 °C. The detailed uniform design optimization methodology was given in previous work 1 . Uniform test design is a computer simulation accomplished in many industrial applications to produce a model presenting the real process performance. Thus, uniform test design U 10 (10 3 ) experiments of 10-run for 10-level of 3 factors for each flotation temperature was done in a randomized order to obtain the optimum flotation reagent system giving the best separation efficiency. The dosage range of each factor was identified through the preliminary tests to be 0.4-2.2, 0.1-1.0, and 0.3-1.2 kg/t for the depressant, activator, and collector respectively. The interactions between the factors in terms of separation efficiency were elucidated using response surface modeling.
Modeling and simulation. GUI of MATLAB was conducted to simulate the optimum flotation reagent system at various flotation temperatures based on the estimated optimum data using a uniform test design. The interface content design of the flotation reagent system calculation program is shown in Fig. 2.
The layout shows the flotation temperature as an input value, while the values of depressant dosage, activator dosage, collector dosage, concentrate yield, grade, and recovery are output after clicking the calculation button which is the execution command button.
The calculation procedure in Code View was as follows:  Furthermore, according to the above procedures, at the site, the flotation temperature of the flotation circuit can be input in the range from 20 to 45 °C, and the flotation reagent system can be obtained to guide the on-site production practice.

Results and discussions
Uniform design factors interaction and optimization. The complexity of the flotation process closedcircuit required utilizing a uniform test design to examine the interaction between the independent factors and to find the optimum separation efficiency. Figure 3 shows the response surface modeling of each two factors affecting the separation efficiency while maintaining the other factors at their optimum values at each flotation temperature. CaO is essential for activating quartz particles to accommodate collector adsorption selectively while depressing the iron oxide particles by starch. For instance, at 20 °C, the change in the CaO or collector dosages insignificantly impacted the separation efficiency, especially at lower starch dosage, but increasing the starch dosage required higher dosages of CaO and collector to maintain the high value of separation efficiency. Whereas for the interaction between CaO and collector, the separation efficiency is minimum when one of them is at the highest level and the other at the lowest level, slightly increases when both are at the lowest levels and is maximum when both are at the highest levels. Increasing the flotation temperature enhances the significant effects of these factors. For example, at 45 °C, to accomplish the maximum separation efficiency, the CaO should increase to the maximum value of 1 kg/t with increasing the starch dosage to 1.8 kg/t at an intermediate collector dosage of about 0.7 kg/t. From Fig. 3, the optimum reagent systems giving the maximum separation efficiency (maximum grade and recovery) at each flotation temperature were identified (Table 1). In this study and a previous study 1 at 30 °C, and at the optimum conditions of 1.6 kg/t depressant dosage, 1.0 kg/t activator dosage, and 0.8 kg/t collector dosage, the maximum Fe grade and recovery of 68.90% and 92.62% respectively were accomplished. In case, the uniform test design is not employed, the typical plant conditions are 1.2 kg/t, 0.5 kg/t, and 0.8 kg/t depressant, activator, and collector respectively; and the plant produces 62.49% concentrate yield containing 68.28% total iron with 89.07% iron recovery 1 .

Mathematical models estimation. Given the optimum reagent systems of closed-flotation circuits at
different temperature values shown in Table 1, the mathematical models were generated for predicting the concentrate yield, grade, and recovery as well as the optimum flotation reagent systems at different flotation temperatures.
The fitting curves of the mathematical model of the optimum flotation reagent system (i.e., depressant, activator, and collector dosages) versus the flotation temperature are shown in Fig. 4. It indicates that the increase in flotation temperature from 20 to 45 °C increased the optimum values of both depressant and activator dosages from 1.0 kg/t and 0.7 kg/t to 1.8 kg/t and 1.0 kg/t respectively, but reduced the optimum collector dosage from 1.2 kg/t to 0.7 kg/t.
At the optimum flotation reagent systems at different flotation temperature conditions, the fitting curves of the mathematical model of concentrate yield, total iron grade, and total iron recovery versus flotation temperature are shown in Fig. 5.  According to the generated models, varying the flotation temperature at the site required adjusting the flotation reagent system with predicting the optimum performance. Thus, the flotation temperature in the flotation circuit can be input in the range from 20 to 45 °C, and the flotation reagent system can be obtained to guide the on-site production practice by giving the expected flotation indexes as shown in Fig. 6.  Table 2. Also, the flotation results obtained experimentally were compared to the values predicted by the calculation procedure as shown in Table 3, and the relative error values between the test values and the data obtained by the calculation procedure are shown in    www.nature.com/scientificreports/ It can be seen from Table 4 that the relative error values are small, and the prediction effect of the calculation program is good. Besides, the adjustment of the reagent system can be quickly guided according to the results of the calculation procedure, and the closed-circuit flotation indexes can be obtained.

Conclusion
In this study, MATLAB-based GUI was utilized to simulate the effect of flotation temperature on the optimum flotation reagent system of the iron oxide reverse flotation circuit. The optimum data were generated first by conducting uniform test design which is a computer simulation utilized in many industrial applications to produce models for presenting real process performance. The increase in flotation temperature from 20 to 45 °C increased the optimum values of both depressant and activator dosages from 1.0 kg/t and 0.7 kg/t to 1.8 kg/t and 1.0 kg/t respectively but reduced the optimum collector dosage from 1.2 to 0.7 kg/t. At these optimum conditions, the product yield, Fe grade, and Fe recovery were in the range of 62.27-63.87%, 67.21-68.90%, and 88.83-92.62% respectively. Then MATLAB-based GUI generated the mathematical models of the optimum reagent systems at different temperatures, so that by varying the flotation temperature at the site, the flotation reagent system can be effectively adjusted with predicting the optimum performance. In summary, GUI is extremely endorsed by mineral processing researchers, and it can also guide the addition of on-site flotation reagents at different flotation temperatures in addition to predicting the flotation indexes such as concentrate yield, grade, and recovery.

Data availability
The datasets generated and/or analyzed during the current study are available for the corresponding author upon reasonable request.